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Abstract 

We propose solar-shaped thermal concentrators designed with orthoradial layers and ob¬ 
tained in practice through the homogenization of an ideal thermal concentrator. Considering 
the spectral regime of the heat equation, we quantitatively evaluate at different pulsations 
the effectiveness of the homogenized concentrators by comparing the thermal flux existing 
in an ideal concentrator and the thermal flux in an homogenized concentrator. Dependence 
on the pulsation is shown to be negligible and plotting the effectiveness of the homogenized 
concentrators as a function of the number of orthoradial layers N, we determine the number 
of layers needed to achieve a certain effectiveness. Significantly high numbers A(ranging 
from a hundred to tens of thousands layers) are found highlighting the fact that achieving 
high effectiveness demands a high level of engineering of the homogenized concentrator. 


1 Introduction 

Following the proposal of Smith and Pendry’s group to design electromagnetic cloaks and con¬ 
centrators via spatially varying anisotropic media deduced from a geometric transform [U, it 
has been shown that thermal cloaks and concentrators could be envisaged using the same tech¬ 
niques [2]. Following these results, recent studies have been conducted to prove experimentally 
the feasibility of a thermal cloak [111!. Additional studies on the control of heat diffusion have 
also shown that considering crossover from thermal diffusive effects to phonon scattering, one 
can design devices for heat management such as thermocrystals [5] and superlattices [6]. 

Regarding transformation optics applied to thermodynamics, an important issue is to en¬ 
gineer metamaterials so as to approach the best we can the required material parameters. A 
homogenization path toward multi-layered cylindrical thermal cloaks with high effectiveness has 
been detailed in recent work [7|. In the present letter, we would like to present a similar ho¬ 
mogenization path toward thermal concentrators with solar-shaped mantle i.e. with a structure 
periodic in the azimuthal direction [8]. We stress that while advances have been made in the 
fabrication and characterization of layered thermal cloak mm, one cannot achieve a thermal 
concentrator with an alternation of concentric layers, as we shall see in the sequel. On the other 
hand, the solar-shaped mantle structure that we here discuss provides an exact solution to ap¬ 
proach the ideal anisotropic concentrator. Ideal and homogenized concentrators are studied in 
the spectral regime and it is showed that the performances of the homogenized concentrators 
do not depend on the pulsation. Furthermore, numerical results allow us to extrapolate on the 
number of layers N that one should implement to achieve a certain effectiveness. 


1 


2 Space transformation 


We consider the radii Ri, R 2 and R 3 such that Ri < R 2 < R 3 and the two-dimensional 
transformation {r', 6 ') = if{r),9) with / defined as: 

f f{r) = ^r, for r E [0,/?2], 

\f(r) = ar + (3, ioi r G [R 2 , R 3 ], 

with a = > 0 /3 = R3^E^ < 0- This transformation maps the region 0 < r < R 2 

onto 0 < r' < Ri (compression of thermal space) and the region R 2 ^ r < R 3 onto Ri < r' < R 3 
(extension of thermal space). 



Figure 1: (a) Normalized conductivity (solid lines) and normalized product pc parameters 
(dashed lines) used for an ideal concentrator. The inner boundary of the concentrator is Ri = 
0.15 mm and the outer boundary is R 3 = 4i?i = 0.6 mm while we have R 2 = 2i?i = 0.3 mm. 
Notice that the behaviour of k'qq is obtained by simply calculating k'qq = (b) Map of the 

temperature in a homogeneous medium and (c) with a concentrator, both at w = 0 rad.s”^. 


The heat equation resulting from this transformation can be written 

BT 

div (k'VT) = (2) 

with T the temperature and the transformed conductivity k' and transformed product of heat 
capacity and density p'c' written as 


^'= kR(0) 


pc'{r) = pc 




where R(0) is the rotation matrix of angle 0 and we have 

/ r df , 1 

^rr = = —• 

/(r) dr k'„ 


(3a) 

(3b) 


(4) 
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Writing r' = /(r), we can derive and k'qq as a function of either r or r'. Considering that 

"^00 = 4 ;’ 

r- € [0, R 2 ] 




for 


r' G [0, i?i] 


and 


r + P 

r' -13 


Rl 


P'c' = pC-^2 


j for r € [i22,^3] 

for r' G [-Ri,i? 3 ] 

r G [0, R2] 


j j 


1 r 


P C pc 2 a 
7- _l_ id 


1 r'-/3 


= 2 / 

r' 


r'G [0,i?i] 
for r G [^ 2 ,-^ 3 ] 

for r'G 


Introducing the Fourier transform of the temperature 

/ OO 

r(x, 

-OO 


(5) 


( 6 ) 


(7) 


where x is the space variable, uj the angular frequency and t the time variable, one can write 
the time-harmonic heat equation as 


div (kS/u) + jojpcu = 0. (8) 

We study this ideal concentrator behaviour for several angular frequencies u by solving the time- 
harmonic heat equation ([5]) using the commercial hnite element software COMSOL Multiphysics 
®. We consider a square computational domain of side length L = 5 x 10“^ m. The medium 
surrounding the concentrator is made of zinc and its thermal properties are k = 121 W.m“^.K“^ 
and pc = 2.75 x 10® J.m“^.K“^. The top temperature is set to Tj = 230°C and the bottom 
temperature is set to Tq = 30°C. Neumann boundary conditions = 0 (perfect insulating 
conditions) are set on the two other sides of our computational domain. 

The behaviour of and p'd for 0 < r' < i ?3 is represented on Fig. UKa). One should notice 
that singularities do not appear in the ideal concentrator inner structure as the terms h'qq 
and p'c' do not vanish or go to infinity in the r' G [i?i, i? 3 ] as it would happen with an ideal 
invisibility cloak. Moreover, the object in the inner region of the concentrator is not chosen at 
will but depends on the surrounding media: its normalized radial thermal conductivity is 
equal to 1 (thus the thermal conductivity of the object is k) and its product p'd is proportional 
to pc by a factor Both the inner region 0 < r' < i?i and the annulus region Ri < r' < R 2 
need to be filled with the according materials for the concentrator to be operational. In the 
end, a structure with such heterogeneous anisotropic characteristics will behave as a thermal 
concentrator as it will concentrate the temperature gradient that would exist on the region 
0 < r < R 2 without the concentrator, to the smaller region 0 < r' < i2i, as it is presented in 
panels (b) and (c) of Fig. [TJ The temperature gradient is actually enhanced by a factor The 
profile of the temperature along the diffusion direction at pulsations a; = 0 rad.s“^ and uj = 100 
rad.s”^ is plotted in Fig. [2j We can clearly see that the gradient of temperature existing on the 
segment [—R 2 ,R 2 ] is constrained on a the smaller segment [—at w = 0 rad.s“^. Fig. [3^ 
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(a) Profile of temperature (b) Profile of temperature 




— Homogeneous medium — Concentrator 

Figure 2: Profile of temperature along the diffusion direction (x = 0 and y E [—2.5 x 10~^, 2.5 x 
10“^]) at (a) a; = 0 rad.s“^ and (b) cu = 100 rad.s“^. The inner radius of the concentrator is 
Ri = 0.15 mm while the outer radius is R 3 = 0.6 mm. Vertical black dashed lines denote the 
boundaries of the inner and outer circles. Vertical green dashed lines denote the boundaries of 
the initial circle of radius R 2 in the virtual space. 


and Fig. [ 8)3 show the variation of the thermal flux by surface unit along the diffusion direction 
at a; = 0 rad.s”^ and cj = 100 rad.s”^ respectively, where the thermal flux by surface unit is 
defined by 

<^ = -(k-VT)-u, (9) 

where k is either a scalar or a matrix depending on the domain and u is the direction vector of 
heat diffusion. With R 2 = 0.3 mm and i?i = 0.15 mm, the thermal gradient is enhanced by a 
factor ^ in the region 0 < r' < i?i. Thus, if we choose a smaller value of Ri, one can achieve 
higher enhancement of the thermal gradient in the inner region of the concentrator as shown in 
Fig. [ 3 ^ and Fig. [ 3 ) 3 . 

Now to obtain such a structure that mimics the thermal concentrator, we design a solar¬ 
shaped mantle in the next section. 


3 Homogenized solar-shaped mantle 


We consider a periodic alternation of layers of respective conductivities ki and K 2 and of re¬ 
spective heat capacities pici and P 2 C 2 (see Fig. HI) so that the overall conductivities and heat 
capacities can be written 


Ke = Ki(x)lffjlj 



PeCe = /3lCi(x)l[o 1] 



+ P2C2{x)t^i^ 



( 10 a) 

( 10 b) 


where x = (xi,X 2 ) is the position vector. 


Arctan 


is the angular coordinate and 1 / 
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— Homogeneous medium — Ri = 0.15 mm — i?i 


(b) Profile of thermal flux 

xlO^ at a; = 100 rad.s“^ 

16 r 



Figure 3: Profile of thermal gradient along the diffusion direction (x = 0 and y E [—2.5 x 
10“^, 2.5 X 10“^]) with i?i = 0.15, 0.1, 0.03 and 0.01 mm leading to a great enhancement of 
the thermal gradient in the region 0 < r' < at (a) a; = 0 rad.s“^ and (b) oj = 100 rad.s“^. 


(Kl.PlCl) (Kz.PzCz) pc''"™) 




Figure 4: Schematic view of the solar-shaped homogenized concentrator. Each element is either 
filled with materials (ki,/9iCi) or {k 2 ,P 2 C 2 ) and those elements are periodically laid out with e 
periodicity. As the number of layers N goes to infinity (or equivalently, e —>■ 0), this structure 
converges to an effective medium of thermal parameters (pc)^"™'. 


is the indicator function of the set I. The function is periodic on the set Y = [0,1] and e 
represents its periodicity (clearly, the thinner the layers in the cloak, the smaller the positive 
parameter e). In this set of orthoradial layers, the time-harmonic heat equation is written 

div (KsVUe) -h jupCeUe = 0 (11) 


where is the Fourier transform of the temperature and oj is the angular frequency of a periodic 
heat source. Following the work from G. Nguetseng and G. Allaire on two-scale convergence 
mm, it is possible to show that when e —>■ 0, this alternation of materials behaves like an 
anisotropic inhomogeneous medium of conductivity 





? ) 


( 12 ) 
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and heat capacity 


1 


{pc)^°^{x) = / pocody = -(pici + P2C2) 


(13) 


Jy 2 

where kq and pqCq are the respective two-scale limits of Ks and peCe when e —>■ 0 and are written 


Ko{x,y) = i](y) + K 2 {x)t^i ^{y) (14a) 

poCo{x,y) = piCi(x)l[o i](y) + p2C2{x)l^i^^^{y) (14b) 

The homogenized time-harmonic heat equation is then written 

div (^A’^°^{x)Vu{x)'^ + Mp{x)c{x)f°^u{x) = 0 (15) 


where u is the limit of when e ^ 0. Then, considering ()14ap and (|14bp . we need to solve the 
system 

r = k ’ 

\ (pc)'^'’™ = pcdet(J) 

+ K 2 ) = Krr 

(16) 




2Kin2 _ 
K,l-\-K,2 


, i(pci-b PC 2 ) = pcdet(J) 


for our set of materials to work as the thermal concentrator. Solving for and K 2 provides us 
with the following solution 




ACl 





(17) 


The equation {pc)^°‘^ = pcdet(J) is not deterministic on the expression of pici and p 2 C 2 - We 
can choose that 


Pi Cl = P 2 C 2 = pcdet(J) = pc 


I r' -(5 


(18) 


r 

The behaviour of ki and K 2 is plotted on Fig. [5^. It is interesting to notice that if one tries to 
obtain a thermal concentrator by homogenizing a set of concentric materials, the solutions of 
m are imaginary numbers. The homogenized matrix with concentric layers is given by 


A 


hom i 


[X = 


= I 


0 ^ 
fyKodyy 


and the system we must resolve becomes 


j^hom _ 

^pc)honi ^ pcdet(J) 


2K-IK2 _ I 

ki+k, 2 rr 


+ K 2 ) = k'qq 

, i(pci-b PC 2 ) = pcdet(J) 


(19) 


( 20 ) 
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Figure 5: (a) Behaviour of ki and K 2 for r' € [0;i?3]. (b) Horizontal configuration of the 

homogenized concentrator with 10 layers: layers are implemented starting from angle 0 = 0. (c) 
Vertical configuration of the homogenized concentrator with 10 layers: layers are implemented 
starting from angle 9 = ^. 


This system gives the imaginary solutions 



Such imaginary conductivities cannot hold as physical conductivities. Thus, it is not possible 
to mimic a thermal concentrator with concentric layers. 

4 Numerical results 

We want to evaluate the ability of one homogenized concentrator to concentrate the thermal 
gradient inside the region 0 < r' < compared to the ideal concentrator resulting from the 
space transformation presented in ([1]). We introduce the effectiveness criterion defined as the 
ratio between the normalized surface integral of the thermal flux in the region 0 < r' < i2i with 
the homogenized concentrator and the normalized surface integral of the thermal flux in the 
homogeneous medium in the region 0 < r < 

—^ ffn ^{homogenized concentrator)dS 

Thermal flux ratio = ^^^--- (22) 

JJr2 TTlcdiUTTljdS 

where is defined in Q. As stated before, the ratio between thermal flux with the ideal 
concentrator in the region 0 < r' < and the thermal flux in the homogeneous medium in 
the region 0 < r < i ?2 is equal to ■^. Therefore, the better the homogenized concentrator, the 
closer the Thermal flux ratio to 
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Thus we performed simulations for homogenized concentrators ranging from N = 3 to N = 
100 orthoradial layers and calculated the Thermal flux ratio obtained in the region 0 < r' < 
at a; = 0 rad.s“^ and u = 100 rad.s“^. We distinguished two configurations for the homogenized 
concentrators: a horizontal configuration and vertical configuration as represented in Fig. [5)3 
and [5):. Horizontal configuration corresponds to the configuration where orthoradial layers 
are implemented from angle 0 = 0 whereas vertical configuration denotes orthoradial layers 
implemented from angle 9 = ^. Depending on the symmetry of the homogenized concentrator, 
heat fluxes won’t experience the same structure. Results are presented in Fig. [6l As we can see, 


(a) Horizontal configuration (b) Vertical configuration 




(c) Horizontal configuration 



Number of orthoradial layers 


(d) Vertical conhguration 



Number of orthoradial layers 


Figure 6: Convergence of the thermal gradient obtained in the region 0<r^<Ri = 0.15 mm 
with homogenized concentrators of orthoradial layers ranging from N = 3 to N = 100 at a; = 0 
rad.s”^ in (a) horizontal configuration and (b) vertical configuration and at cu = 100 rad.s”^ in 
(c) horizontal configuration and (d) vertical configuration. 






























both horizontal and vertical configurations of the homogenized concentrators give oscillations 
of the thermal flux ratio. These oscillations appear both at a; = 0 rad.s”^ and oo = 100 rad.s”^ 
and converge towards the same limit This highlights the fact that the convergence of the 
Thermal flux ratio is independent of the angular frequency io (other values of cu would give 
the similar oscillations converging towards the same limit). This result is to be expected since 
the angular frequency uj does not appear in the homogenized parameters as we can see from 
equations (fT5]) . (fT71) and (fTSl) . Now, independent of lo, even values of N in horizontal and vertical 
configuration give significantly the same value of Thermal flux ratio for a given number N of 
layers whereas odd values of N lead to different Thermal flux ratio in horizontal and vertical 
configurations. This result is explained by the fact that, for even values of N, the homogenized 
concentrator is perfectly symmetric with respect to the heat diffusion direction, in horizontal 
as in vertical configuration. As a comparison, odd values of N do not ensure the symmetry 
of the homogenized concentrator. These results illustrate the fact that the effectiveness of 
homogenized concentrators depends substantially on the layout of the orthoradial layers. Such 
considerations of symmetry do not appear with the homogenization of cloaks since these involve 
only concentric layers meaning perfect symmetry whatever the direction of the heat diffusion. 
We see however that the difference of oscillations between even and odd values of N tend to 
zero as N increases as one should expect. Now to understand further the behaviour of the 

(a) Heat Flux convergence (b) Heat Flux convergence 




Figure 7: Heat flux convergence as a function of log(N) for the horizontal configuration at (a) 
a; = 0 rad.s“^ and (b) (a) io = 100 rad.s“^. Linear fitting has been applied to all values of N, 
and even and odd values of N separately. 


homogenized concentrators, we define the quantity 


Heat Flux convergence = 
1 


ttRI 


// ^{homogenized concentrator)dS — // ^{perfect concentrator)dS 
JjRi JjRi 


(23) 
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This quantity is going to zero as e —)• 0 (i.e N going to infinity). We can use the result [13] 


||tt — rie||L2 ~ Y y 1^ “ rtepdx < C(n)e (24) 

where rtg denotes the temperature in the computing domain O for a Wlayer concentrator and 
u is the limit of when e —>■ 0. Plotting the decimal logarithm of the Heat Flux convergence 
as a function of log(A^), we can verify if the behaviour of Heat Flux convergence as N goes 
to infinity is consistent with (|24ll for both a; = 0 rad.s“^ and ui = 100 rad.s“^. Results are 
illustrated in Fig. [T] 

Linear fitting has been applied to all values of N, and even and odd values of N separately. 
We find that 

log{Heat Flux convergence) ks —log(A^) + C (25) 

for the three sets of points in Fig. [7[i and Fig. (TJr, confirming the behaviour expected from (1241) . 
As seen above, the homogenization convergence is again independent of angular frequency oj as 
one should expect. Using the linear fit for all N values (blue curve), we can extrapolate to find 
the required number N of layers to achieve a desired value of Heat Flux convergence. 

The relative sensibility is defined by the ratio of the Heat Flux convergence and the 
normalized heat flux ffj^^ ^{perfect concentrator)dS. The number N of layers as a function 

of relative sensibility for both w = 0 rad.s“^ and uj = 100 rad.s“^ is plotted on Fig. [8] in the 
case of horizontal configuration. We see for instance that to reach a relative sensitivity of 10“^, 
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Figure 8: Number of layers as a function of relative sensitivity at cu = 0 rad.s ® (blue) and 
LV = 100 rad.s“® (red) in the horizontal configuration. 

one would need to implement a homogenized concentrator with nearly 2000 layers. A relative 
sensitivity of 10“^® would need an even higher number of layers, close to 20000. Therefore, we 
can determine the order of magnitude of the number of layers that is needed to reach a certain 
effectiveness. However, the order of magnitude of the number of layers on Fig. [8] highlights the 
fact that achieving high effectiveness demands a high level of engineering of the homogenized 
concentrator. 
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5 Conclusion 


While concentric multilayers do not allow to reach perfect heat concentrators, it was shown 
how the homogenization of orthoradial multilayered structures allows one to achieve an ideal 
anisotropic heat concentrator. The complete design of such structures was given and we cal¬ 
culated the temperature and flux variations within the thermal device at different pulsations. 
Furthermore, the quantitative effectiveness of the homogenized concentrators was numerically 
calculated versus the number of orthoradial layers and it was shown that the performances of the 
multi-layered concentrators are independent of the pulsation. Extrapolation on the numerical 
results allows one to adapt the complexity of the concentrator once the required performances 
are chosen and fairly high numbers of layers are found. This study emphasizes the fact that en¬ 
gineering devices for the heat management requires optimal designs but paves the way towards 
highly effective thermal metamaterials. 
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